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Abstract: Primordial cosmological hypermagnetic fields polarize the early Universe 
plasma prior to the electroweak phase transition (EWPT). As a result of the long range C ^ 
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plified, opening a new perturbative way of seeding the primordial Maxwellian magnetic ^ — >^ 
field at EWPT. 
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The electroweak phase transition (EWPT) has long been considered as a playing an 
important role in the generation of primordial magnetic fields ||l], ||, ^. Here we suggest 
that the seed fields accounting for the observed intergalactic magnetic fields arise, say, 
as tiny fluctuations ~ 10~^^ Gauss, associated to some early phase transition prior to 
the EWPT. The interplay between the resulting polarization effects of the early Universe 
plasma and the long range parity violating gauge interaction present in the Standard Model 
subsequently amplifies the seed field till an epoch close to the EWPT time, after which its 
evolution is described by standard perturbative physics, free from uncertainties from the 
EWPT epoch. 

It is well-known that, at the high-temperature symmetric phase of the Standard Model 
(SM) all gauge bosons acquire a "magnetic" mass gap ~ g^T, except for the Abelian gauge 
field associated to weak hypercharge. Collective plasma effects also modify the dispersion 
relations of the quarks and leptons at finite temperature giving them effective "chirally- 
invariant" masses determined by their gauge group quadratic Casimirs Q. 

While long-range magnetic-like fields associated to non-Abelian gauge forces do not 
exist, the Abelian hypercharge magnetic fields are never screened and can survive in the 
plasma for infinitely long times. The nonconservation of the lepton charge due to the 
Abelian anomaly associated with the hypercharge field has been suggested to play a 
role in generating the observed the baryon asymmetry of the Universe . 

Here we stress the importance of polarization effects in the early Universe plasma 
associated to the parity violating weak hypercharge interactions well before the EWPT. 
We show how, in the presence of a lepton asymmetry expected to arise in the leptogenesis 
scenario [^], a pre-existing field B^^ , for simplicity assumed to be a large-scale field, will 
undergo spectacular growth prior to the EWPT. Below the critical EWPT temperature 
the hypermagnetic fields are converted into standard Maxwellian magnetic fields. Thus the 
amplification mechanism can seed the ordinary magnetic field, helping to account for the 
observed magnitude of intergalactic magnetic fields Q . 

Consider the equations of motion for the hypercharge l^-field in the hot plasma and in 
the presence of a pre-existing large-scale hypermagnetic field B^, regular on scales smaller 
than the horizon size at Tq > T Tewpt- We assume that this primordial field has a 
small amplitude g' B^^ <C T^. For simplicity, we neglect the Abelian anomaly and assume 
flat Minkowski space. The U{1)y interaction Lagrangian in the SM with the one Higgs 
doublet (f'^ = (/)('^)) is given as. 



1- 1- 4- 2- 



+ 



a 

+i^—-[^+D^^-{D^^+)^]. (1) 

For the assumed external seed hypermagnetic field = V x Y*^^-* = (0,0, i?^), this leads 
to the following Dirac equations for massless charged leptons = e, /i, r ), neutrinos and 



quarks {Ui = n, c, t, Di = d, s, b) 



_/{a)(^')y(o)l ^(a) 

where a = Il,Ir,v^,Ul^ Dl,Ur, Dr and f^"'\g') denote the corresponding SM couplings: 
fRig) = ~9' for right-handed charged leptons, fiig') = —g'/'^ for left-handed charged 
leptons and neutrinos, f^'^\g') = g' /Q and f^^\g') = 251' /3, f^^\g') = —g'/'i for left- 
handed quarks and right-handed quarks respectively. 

The resulting Landau spectrum in the mean hypercharge magnetic field takes the form 



e{pz,n, A) = y Pz+ \ fL,R{g') \ B(^{2n -|- 1 =F A), where the upper sign applies to particles 
and the lower one to antiparticles, irrespective of their chiralities, A = ±1. Within the 
small-field approximation g Bq <^ we have 

eip,X)=p+\fL,R{9')\B^X/2p (2) 



with p = yp^+p\ and p\ =| fL.R.{g) I BQ^{2n + 1). Note that such paramagnetic term 
in Eq. (|2|) coming from the spin of fermions is absent for Higgs bosons. 

The above spectrum leads to the unique definition of the equilibrium distribution func- 
tions in the multi-particle approach. Up to the Abelian anomaly lepton (and quark) num- 
bers are conserved through the continuity equations djjj^^ / dx^ = Qioi a = l^, £r, ui, U l,/?, Dl^r. 
Therefore we may define chemical potentials ^^""^ for the corresponding equilibrium distri- 
bution functions. From Eq. (|2|) we easily find the equilibrium one-particle density matrix 
for massless particles 

^A'A exp[(e(p,,n,A)^//„)/r] + l 
which is approximated in the quasi-classical limit, by 

/yf ^ ^/i"''^ W + ''-^s^r^'ipl (3) 

with the lower sign for chemical potentials corresponding to antiparticles. Here = 
[g(pT/^a)/7' _j_ g^pg ^Y^Q Fermi distributions of particles (antiparticles) which for iXa/T <C 1 
correspond to the Lorentz-invariant densities n^"'"\ 
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(27r)3exp(b^^J/T) + l 

(4) 
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Here Ueq = 3C(3)T^/47r^ is the equilibrium lepton density when = 0; ("(3) ~ 1.202 is 
the Riemann function value. The second term in Eq. (^) includes Pauli matrices multiplied 
by the mean spin vector, 

St'\p) = = boS^^)(p), (5) 
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which is the equihbrium spin distribution function. Here = /Bq is the unit vec- 
tor along the mean hypermagnetic field. Notice that equilibrium spin distribution values 
iSq"'"^ (p) define the densities of fermions populating the main Landau level, n = 0: 



(2vr) 

Ug')\B^Tln2 

4^2 



1±— ^ + 0( 

2T\n2 \\T 



Note also that the subindex A in the density matrix corresponds to the spin projection on 
the hypermagnetic field, {(Tz)x'x = -^i^a'a- Together with chirality 75^6^^ = i^e^L is a 
good quantum number since [75, S^] = 0. 

Note that the density asymmetry for massless fermions coming from Eq. , n*^'^^ — 
fii"') = ^^t2/6, is one half that of the Higgs bosons (^aT^/S). Taking also into account the 
equilibrium conditions for the chemical potentials of the plasma components (using plasma 
neutrality < Q >= 0, < Q3 >= fiw = 00) and normalizing the hypercharge Y on T^/G 
where Y directly reads from the Lagrangian Eq. (||) as ^, 

Y = - ^iHuiL + /ii) - 2 ^ /i'^ + N[fi^L + /UdL + 4/U„i? - 2fidR] + 2(^+ + ^0) = 
I I 



2(Q-Q3) = 2 



-2^^'i + 6^„L + 14^o 



0, (7) 



one finds the chemical potential of the neutral Higgs boson /Uq , 

/^O = j • (8j 

Using the sphaleron equilibrium condition valid above EWPT, 

Y,l^i = -^f^uL (9) 
I 

and the expressions for the baryon and lepton numbers in terms of chemical potentials 
(given in Eq. (8) in Ref. ^) one can see that /Uq vanishes only in the exceptional case 
B = L = (see, e.g. in [0]). 

Note that in the first line of Eq. (^) we have used /i+ = — /i- = ^o- This follows from 
the equilibrium relation fiw = fi- + fiQ (implied by conversions W~ <j)^) and the 

neutrality condition for the isospin component Q3 ~ fiw = 0. We also took into account 
the common color factor of quarks and assumed that quark asymmetries are degenerate in 
flavor for all left-handed up-quark fields and all right-handed down-quark fields respectively 
(note that we do not assume such degeneracy for charged leptons and neutrinos). 

One can easily write down Maxwell-like equations of motion for the Abelian Y^-field in 
the equilibrium plasma. These involve thermal averages of the various four-current fermion 
densities, e. g. for leptons < iaj^i^a given as differences of Wigner distribution functions 



""^The same result follows from the definition of hypercharge Y in Eq. (8) in paper 



for particles and antiparticles, 6f[^\ (p, x, t) = f[l'l (p, x, t) — f['^\ (p, x, t). Substituting the 

A A A A A A 

general distribution function, 

/S(p,x,i) = ^/('^)(p,x,t) + ^-^St\p,^,t), (10) 

and summing over spin variables A, A' one can recast the Poisson and Maxwell- like equa- 
tions as 



V • By = 0, V • Ey = 47r t) + Jo^5(x, t)] 



dt 

dEy 



-- -V X Ey, 

+ V X By = 47r [J^(x, t) + J^x, t)] 



These describe the equilibrium plasma at T ^ Tew pt and differ from the familiar Maxwell 
equations by the presence of the pseudovector current associated to the parity violating 



1^ interactions. 



The vector current is a sum = jj^ + 3A^ J^^-^^ + J^^p^^ where each term is given 
in terms of differences of current asymmetries 



xAa) _ -(a) _ -(a) 



(27r)3 p 

in the form following from the interaction Lagrangian (||), 



Jj;(x, t) = -^ 25jff (x, t) + (x, t) + (x, t) 



9_ 
12 



(x, t) + 6jr> (x, t) ) + A6j}r^> (x, t) - 26 j^' (x, t) 



ML), 



(uR), 



AdR) I 



involving partial current asymmetries given as 



(12) 



(13) 



where n^""^ = n('*)(x, t) = jjf^ (x, t)u'^(x, t) are the densities defined in Eq. (|^ which 
coincide with those calculated in the medium rest frame, V = 0, 7 = (1 — y^)^^/^ = 
1. Summing the individual contributions in Eq. ([T2|) one obtains the hypercharge vector 



current density as 



Y 



Uc(3)j 



4T 



tY 



E 



(14) 



As seen from Eq. (0) hypercharge neutrality of the plasma < Y >= implies the 
vanishing the component Jq = at large scales exceeding the Debye radius, r ^ rn, 
resulting in the relation between chemical potentials Eq. (|8|). 

We now turn to the axial current. Using Eq. ([lOD one can express the axial vector 



current JY^{'x., t) in terms of differences of spin distribution functions entering the axial four- 



vectors MJT^ (p, x, t) = <^ (p, x, t) - <^ (p, x, t) with A),'> (p, x, t) = {p- S^'^) (p, x, t)); (p • 
S(-)(p,x,t))pM 2. 

The total pseudovector current = J2i "^i^^sl^j + 3-^<^^5 is given in terms of partial 
contributions of leptons and quarks: 



J, 



£At5 



4 
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■54"^)(p,x,t) 
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■M(r^)(p,x,t) 



3 ./ p(27r)3 M ^t-' '"^ Q J p(27r)3 
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d'^p 



d^p 
p(27r)' 



■M(f^)(p,x,t) + 



M|f^)(p,x,i) 



(15) 



Substituting the equilibrium spin distribution Eq. (||) into Eq. ([T^) we immediately get 
Jo5 = in the rest frame of the medium as a whole, V = 0, 7 = 1. Thus, under the 
hypercharge neutrality condition in Eq. (|^, one has that Jq = Jq^ = 0, so that the 
Poisson equation takes the standard form, V • Ey = 0, in agreement with Ref. ||5|. (Note, 
however, that small scale fluctuations described by non-equilibrium distribution functions 
(5/('')(p,x,t), (5S(")(p,x,t) could lead to 5J^ ^ 0, (5Jo5 0). 

In contrast, each pseudovector 3-current or J^g in Eq. (|l5|) is nonzero even in 
equilibrium, 



eq 



•m- 
J5 . 



+ 4 [Js 



q5Jeg 
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J5 



+ 3 [J5 
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J5 



J5 



(16) 



Substituting the equilibrium spin distributions in Eq. (^) and using again the hypercharge 
neutrality condition and the sphaleron equilibrium condition Eq. (^) we get, after summing 
over leptons J2e quarks, 

/2 



(jr) 



9 



eq 



967r2 



-2 J]]/x^ + 10/z„L + 32^o 



B 



Y 



47 ,0 „y 



15127r2 



(17) 



which leaves then only one independent lepton asymmetry which we take as that of neutri- 
nos, fJ-i = '^ifi'i/iL = /"i^- Notice that we have used here Eqs. (P), (^) for the chemical 
potentials /xq and fiuL- 

^This is the massless limit of the four- vector spin distribution that generalizes the Pauli-Lubanski spin 
vector in the multi-particle approach jl^], A^^'' {p,x,t) — rria x = limm„-.on^a x [pS'^'/ma; S*^"-* + 
(pS'''')p/ma(ep + ma)], where S'"' = S*'^'(p,x,t). Note that the standard Lorenz-invariant normalization 
a^a^ — — 7^ is replaced by A^^^A'^''"'^ = in the massless fermion case. 
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Thus, the system of Magneto-Hydro-Dynamics (MHD) equations for hyperelectromag- 
netic fields in Eq. (11) finally gets the form: 



V By 

dBy _ 

~dr ~ 



0, V • Ey 

-V X Ey, 



47 



dt 



+ VxBy = 4vrJ^ + -^x^Br 
378 vr 



>o 



(18) 



Averaging the total field By = B^ + by(x, t) over random small-scales, < by >= 0, 
we can rewrite Eq. ( p!8| ) as an MHD system for mean hypermagnetic fields By = B^ 
completed by the Ohm equation. In the rest frame V = of the isotropic early Universe 
plasma we are considering the Ohm equation reduces to 



tY 



(19) 



One sees that the MHD equations derived here in the standard Weinberg-Salam model 
qualitatively coincide with what one obtains using the anomaly term in the interaction 
Lagrangian [||, 0]. 

By combining the Ohm law Eq. (|l9| ) and the last Maxwell-like equation in Eq. ( |l8|) 
then using dBy/dt = —V x Ey, we can write the Faraday equation describing the so-called 
a^-dynamo of hypermagnetic field as 

dBy 



dt 



V X aB + ??V^Bi 



(20) 



where r] = (47r(Tcond)~^ is the magnetic diffusion coefficient and we neglect, as usual in 
MHD, the displacement current d^y/dt and use the rest frame condition V = 0. The 
parameter a is the hypermagnetic helicity coefficient given as 



a 



47g ^fi^ 



(21) 



and plays crucial role in the evolution of the hypermagnetic field. We can solve Eq. ( pO|) 
through Fourier harmonics as By(x, t) = f (d^k/ {2tt)'^ By 0^,t)e'^^^ where By{k,t) is ex- 
pressed as 



By{k,t) 



exp 



[a{t)k-i^{t)k^)]dt 



to 



(22) 



For < < a/?7, or correspondingly correlation length scales rj/a < A < oo such field gets 
exponentially amplified, but differently for different scales A. E.g. for the Fourier mode 



k = a/2r] (or A ~ 2ry/a) one gets the maximum amplification 7 = ak — rjk'^ = /Arj ||^, 12| 



By{t) = S^exp 



to 



4rj{t' 



-dt 



or 



By{x) = Bq exp 



32 



dx / ^1.(3; ) 



X 



0.001 



(23) 



where we introduced the new variable x = T/Tewpt and Bq is the assumed initial am- 
plitude of the hypermagnetic field at Tq ^ Tewpt- In the second equality we substituted 
T] = (47rl37T cos^ 0vf)~^) the analogue of the magnetic diffusion coefficent taking into ac- 
count the change from the standard QED fine structure constant (e^ ~ 137^^) to the analo- 
gous hypercharge one, with — > g"^ = e^/ cos^ 6w, where sin^ 6w = 0.23 is the electroweak 
mixing parameter. We have also used the appropriate cosmological time-temperature re- 
lation t = [3.84 X 10^^(r/MeV)~^/\/^ MeV~^ with the number of relativistic degrees of 
freedom g* ~ 100. One sees that, even for small values of the lepton asymmetry one 
obtains a very strong amplification of the seed hypermagnetic field Bq at this scale. 

Unfortunately, the unknown dependence of S,u{x) for x = T /Tewpt > 1 prevents us 
from a reliable numerical estimate for the amplitude By{x). Moreover, since the evolution 
of i^jy depends on nonperturbative physics, it is not easy at this stage to confront it with 



the primordial nucleosynthesis bounds |14]. However, in order to survive against ohmic 



dissipation due to finite conductivity r] = (47r cicond) ^ we should have A > Idiff, where 
^rjln is the diff'usion length. This leads to an upper bound on (,uix) = ^y{T)/T, 



•■diff 



0.001 



x>l, 



(24) 



which explicitly depends on the chosen scale A. Hence for the mode Ai = 2ri/a one gets 
A = 0.23 so that the amplification factor is ~ 32 as seen in Eq. (p3|), while for a larger 
scale, say A2 = 16r]/a, one gets A = 1.83 with a reduced growth factor ~ 7.4 coming from 
the general solution in Eq. (p2|): 



Brit) = Bq exp 



256 



Bq exp 



7.4 -r^ 



dx I S.u{x 



0.001 



(25) 



In any case one can have very strong amplification even for larger scales for which the 
factor ^(A) is bigger. 

In contrast to the mechanism suggested in Refs. and Q ours does not rely on 



the Chern-Simons anomaly term in the SM Lagrangian [ 11 1 . The presence of the anomaly 
acting at the later EWPT epoch, could play an important role in the subsequent evolution of 
the lepton asymmetry produced by the parity violating hypercharge interaction. However, 
by then the asymmetry has already induced the strong Maxwellian magnetic fields which 
no longer convert to leptons, as these carry no anomaly, their evolution being described by 
standard MHD equations. Here we do not study the EWPT conversion of hypercharge field 
to the Maxwellian magnetic field B. However we note that the seed value By ~ 0.3T'^ < 
'^EWPT ~ ^^^'^ Gauss can be easily reached through our Eqs. (23), (25). This provides 
a strongly first order EWPT that, in turn, allows to avoid the sphaleron constraint for 
baryogenesis within the Standard Model [|l^] . It is also important to note that our bound on 



the neutrino asymmetry in Eq. (24) provides the large scale of the mean hypermagnetic field 
A ~ Lq in Eq. (1) of ref. [15|, so that bubble formation during EWPT takes place in the 



background of essentially constant field |15|. The subsequent evolution of the Maxwellian 



field proceeds through the inverse cascade |16| 



Let us now comment on the physical interpretation of the new magnetic heUcity term. 
The original seed field Bq polarizes the fermions and antifermions (including neutrinos) 
propagating along the field in the main Landau level, n = 0. This polarization effect 
causes fermions and antifermions to move in opposite directions with a relative drift velocity 
proportional to the lepton asymmetry. The existence of a basic parity violating hypercharge 
interaction in the SM induces a new term in the hypermagnetic field in Eq. (^) V x aBy 
which winds around the rectilinear pseudovector hypercharge current J5 parallel to By. 
This term amplifies the seed hypermagnetic field according to Eqs. (^) and (|23|). It is 
interesting also to consider in detail the hypermagnetic helicity H = J d^xYBy and how 
it becomes the magnetic one below the EWPT |17]. 

In summary, while we are far from having a complete and fully quantitative picture 
for the origin of intergalactic magnetic fields, we think that the amplification mechanism 
described here could play an important role towards the goal of accounting for the observed 
intergalactic magnetic fields from first principles. 
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